Spin dynamics from Majorana fermions 
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Using the Majorana fermion representation of spin-1/2 local moments, we show how it is possible 
to directly read off the dynamic spin correlation and susceptibility from the one-particle propagator 
of the Majorana fermion. We illustrate our method by applying it to the spin dynamics of a non- 
equilibrium quantum dot, computing the voltage-dependent spin relaxation rate and showing that, 
at weak coupling, the fluctuation-dissipation relation for the spin of a quantum dot is voltage- 
dependent. We confirm the voltage-dependent Curie susceptibility recently found by Parcollet and 
Hooley [Phys. Rev. B 66, 085315 (2002)]. 
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The mathematical difficulties of representing spins in 
many body physics have long been recognized. The 
essence of the problem is that spin operators are non- 
abelian: they do not obey Wick's theorem and an ex- 
pectation value of the product of many spin operators 
cannot be decomposed into products of two-operator ex- 
pectation values, even within a free theory. 

A conventional response to this difficulty is to represent 
spins as bilinears of fermions Q or as bosons . One of 
the disadvantages of these approaches is that the Hilbcrt 
space of the fermions or bosons needs to be restricted by 
the application of constraints 0. [J Q Another difficulty 
is the "vertex problem" , which arises in the context of 
spin dynamics and spin relaxation. Once the spins are 
represented as bilinears, the spin-spin correlation func- 
tions are represented by two-particle Green's functions. 
The calculation of these quantities requires a knowledge 
of both the four- leg vertex and the single-particle Green's 
function. Typically, the vertex is simply neglected, or 
treated in a very approximate fashion. 

An alternative approach is to take advantage of the 
anticommuting properties of Pauli matrices, writing the 
spin operator in terms of Majorana fermions 0, SH IS 



S = --jpx if, 



(1) 



where if = (r] , rj 2 , rf) is a triplet of Majorana fermions 
which satisfy {f? a ,?7 b } = 8 ab . This representation does 
not require the imposition of a constraint: the fact that 
S 2 = 3/4 follows directly from the operator properties of 
the Majorana fermions. In this letter, we show how this 
representation also solves the vertex problem. To demon- 
strate this, we employ an alternative derivation of the 
Majorana spin representation. Consider a spin-1/2 op- 
erator S with dynamics described by a Hamiltonian H. 
Let us now introduce a single Majorana fermion $ which 
lives in a completely different Hilbcrt space, commuting 
with S and H . It follows that $ is a fermionic constant 
of motion, d$/dt = —i[H,$\ = 0: an object of fixed 
magnitude < I >2 = 1/2 which anticommutes with all other 



fermion operators. We may now identify if in Q with 
the operator identity 



if = 2<$>S. 



(2) 



We may confirm that {rj a ,rj b } = S ab using the anticom- 
muting algebra of spin-1/2 operators {S' a ,S' b } = \5 ah . 
Furthermore, using the SU(2) algebra of spins, if x rj = 
2S x S = 2iS, from which Q follows immediately. As 
a last step in the derivation, we note that the indepen- 
dent Majorana operator can also be written in the form 
$ = — 2ir] 1 rfrf, an object that can be verified to com- 
mute with expression 0). (Notice, incidentally, that al- 
though it is true that S = <f>if, this expression is of limited 
use because if and $ are not independent fermions: they 
commute, rather than anticommuting.) 

The important, yet previously unemphasized feature 
brought out by this derivation is that the Majorana 
fermions and the spin operator are proportional to one 
another, if cx S, where the constant of proportionality is 
a constant fermion. From this fact, it follows that 



lv a (h)v b (h) = 2$(t 1 )$(i 2 )^(i 1 )S b (i 2 ) 



S Q (ii)S 6 (t 2 ). 



(3) 



This operator identity enables us to connect the spin 
correlation function to a one-particle Majorana Green's 
function. Inserting commutators or anticommutators 
into ) and taking the expectation value, we find 
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The expectation of a spin anticommutator is a correlation 
function, but its fermionic counterpart represents a re- 
sponse function. Likewise, the expectation value of a spin 
commutator represents a spin response function, but this 
is equal to a fermion correlation function or "Kcldysh" 
Green's function |^. Thus the correlation function of 
the Majorana fermions determines the response function 
of the physical spins, and vice versa. 

We may formalize this relationship, writing 



x"(h,t 2 ) = - A G K (t x ,t 2 ) 



G R (t 1 ,t 2 ) — G A (ti,t2) 



where 



WTitiM) = -{[s a {h),s b {t 2 )]) 
(C) ab (h,t 2 ) = ^<OTi),S 6 (t 2 )}) 



(4) 
(5) 

(6) 
(7) 



are the spin response and correlation functions and 

Gf(h,t 2 ) = -i{[r, a {t 1 ),r } b {t2)]) (8) 
GtdiM) = -t({v a (ti),v b (t2)})0(h-t 2 ) (9) 
Gf(ti,t 2 ) = i{{r) a {t l ),rf>{t 2 )})6{t 2 -t x ) (10) 

are the Keldysh, retarded and advanced Green's func- 
tions of the Majorana fermion. For most purposes, we 
are interested in systems that are in thermal equilibrium, 
or that have reached a non-equilibrium steady state, for 
which the correlation and Green's functions are functions 
only of the time difference t\ — t 2 . In this case, we may 
transform Q into frequency space, writing 



In equilibrium, the function h(ui) = /ib(w) = 
tanh(u/2T) is determined by the fermionic fluctuation- 
dissipation theorem. We recover the conventional bosonic 
fluctuation-dissipation theorem as the inverse of /ie(u;): 



1 



(x")- + M h E (w) 



coth 



2T 



In non-equilibrium steady state conditions, h(u>) must 
be computed from first principles, as a non-equilibrium 
fluctuation-dissipation theorem, but the inverse relation 
between the spin and fermionic fluctuation-dissipation 
functions is preserved. 

We can apply the Kramers-Kronig relation to deter- 
mine the full dynamic susceptibility from l|13f) . as 



x~ + M 



dv 



A(y)h(y) 



(16) 



so that the static transverse susceptibility is given by 



1 



X± = oX~ + (0) 



dv 



A(v)h(v) 
2v 



(17) 



Since (x") + is the Fourier transform of the response 
function %{[S~(t), S + (0)]), using [£-,£+] = -2S Z , we 
can evaluate the ^-component of the magnetization from 



(t = 0) = -i J dwA(uj)h(u 



(18) 



Thus from the fermion propagator one can read off both 
the spin dynamics and the static magnetization. 

To illustrate this method in its simplest form, consider 
a spin-1/2 in a magnetic field in the negative z direction, 
for which H = BS Z . Written in terms of fermions, 



X» = ^G k (lj), 



Gjj(w) -G a {uj) 



(11) 
(12) 



Here, x"( w ) — ImfX-R^)] — I m [x( w + is the imagi- 
nary part of the retarded spin susceptibility. 

It is particularly useful to combine the Majorana 
fermions into a conventional Dirac fermion, writing f' = 
^(r? 1 + i?/ 2 ), for which {/, /T} = 1. The /-fermion is 

directly proportional to the spin raising operator f> = 
y/2$S + , so that S-(i)S+(0) = f(t)P{0). Recasting the 
steady state version of @ in terms of the raising and 
lowering operators, and Fourier transforming the result- 
ing expressions, we obtain 



(xT + H 



-A(u)h(uj) 



(13) 
(14) 



wheieA=£(G R -G A ) 
function and 



ilmG^(w) is the / spectral 



h(u) 



Gr{uj) - Ga(w) 



(15) 



h = -inWB = b fpf - ~ 

The retarded /-Green's function is given by 



(19) 



(20) 



u- B + iS' 

so that A{lo) = 6(oj — B) 1 and from (|13[) it follows that 



(x'r+M 

C-+(oj) = tt5{uj-B), 



,t,>U' - B)\ i xub ( ^) , 



(21) 
(22) 



and from |JT%)|. (S z ) — — itanh^, recovering the Bril- 
louin function. 

The utility of the method is its ability to handle both 
equilibrium and non-equilibrium situations. To illustrate 
this point, consider a spin coupled to two conduction 
seas, according to the Kondo Hamiltonian 



H = 



H\\> — 



e k\C kXa Ck\a 

k\(T 



A, A' 



XX' 



/ J JxX'C^kXa^afiCk'X' (3 1 S 
k,k' 



(23) 
(24) 
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FIG. 1: Leading order Feynman diagram for the Majorana 
fermion self-energy £ a &(u;), which is isotropic in spin indices 
in zero applied field. The self-energy of the /-fermion f' — 
(V 1 + i-V 2 ) is given by E(w). 



where the terms H\\> (A = L, R) describe the electron 
"co-tunneling" between lead A and A' that is mediated 
by spin exchange with the spin S. This model has been 
used to describe the low energy physics of a quantum 
dot. Even when perturbation methods are applied to this 
model, it is difficult to directly extract the spin dynamics. 
The Majorana method permits the spin dynamics to be 
computed perturbatively in the couplings J\y, without 
any approximation to the spin vertex, even when the two 
leads arc at different voltages. 

The Green's functions in zero field are now given by 



Gr,a( w ) 



1 



G r {uj)^ k {uj)G a {uj), 



where T, RA x(uj) are the retarded, advanced, 
Keldysh self-energies of the /-fermion, so that 



h(co) 



T, r (gj) - Ea(w) 



(25) 

(26) 
and 

(27) 



reflects a change in the fluctuation-dissipation theorem. 
The spin-spin correlation function is thus given by 



Writing E r (uj) = E'(w) 
is Lorentzian, where Z 



1 



1 



(28) 



Ea(w) to - Sh(w) 
— iT(oj), then at low frequencies 
Z T 



2 uj 2 + 



r 2 

1 



(29) 



and r 



o 



(I - dZ'/dw)- 1 
ZT(0) is the spin relaxation rate. 

To leading order in the coupling constant, we can com- 
pute the non-equilibrium self-energies by a simple trans- 
formation of the equilibrium self-energies. (The relevant 
Feynman diagram is shown in Fig. In equilibrium, 
the retarded self-energy, computed by an analytic con- 
tinuation of the imaginary time self-energy, is given by 



E fl H=4^(J AV ) 2 T(a 



(30) 



A, A' 



where J\y = pJ\\> , and p is the density of states per 
spin per lead. 



T(uj) = uj 



Y \2mTJ \2itT J uj 



(31) 
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FIG. 2: Voltage-dependent spin-relaxation rate, computed 
to leading quadratic order in the Kondo coupling constant, 
showing the zero field, Korringa component, and the the 
voltage-dependent inter-lead component. For the case cho- 
sen, J\\i = J (A, A' G {R,L}). 



where V'( z ) = d\nT{z)/dz is the digamma function and 
D is the bandwidth of each lead, which has been intro- 
duced as a Gaussian cut-off in the Feynman diagrams. 
The imaginary part of l|31|l is given by 



Imiy = — coth 

v ' 2 



UJ 



2T 



so that the equilibrium Keldysh self-energy is given by 
Ejf (w) = (E fl - E A )h E (u>) = -4^^(J AV ) 2 w. (32) 

A, A' 

The effect of applying a voltage to the leads can be incor- 
porated into the self-energies by noting that a voltage is 
equivalent to a gauge transformation on the conduction 
electron fields, Ck\ a — > C]~\ G e HLxt : , where \i\ is the chemi- 
cal potential shift in the A = L,R lead. In a second-order 
calculation, this gauge transformation can be incorpo- 
rated by making the replacement ui — > uj — (fi\ — py ) in 
each term of the self-energy, so that at finite voltage 

£ fl (u;) = 4^(J AV ) 2 T[^ - (/i A - p y )}, (33) 



A. A' 



Ejf(w) = -4irtuj[J 2 RR + J 2 LL + 2J 2 RL ], 



(34) 



where the voltage dependence cancels out of the second 
expression. Using l|29l) we may immediately read off the 
leading order expression for the voltage-dependent spin 
relaxation rate of a quantum dot: 



4ttJ 2 



RL 



2aT + eycoth 



2TJ 



(35) 



where we have introduced eV = \Il — M_r and a = 

(Jll + Jrr)I 2 Jrl- ( See Fi S- H) We recognize the 
V = limit of l|35|) as the Korringa relaxation rate of 
a single spin [13|. At finite V, the second term gives the 
voltage-dependent spin relaxation rate induced by the 
coupling between leads: this term is linear in tempera- 
ture for eV <C T, but linear in voltage for eV ^> T. An 
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FIG. 3: Voltage-dependent spin correlator C(uS) calculated 
for a sequence of voltages V = 0, 0.3, 1 at a temperature T = 
0.05. Here a = 1 and J = 0.2. Inset, fluctuation dissipation 
function h(u>) for the same sequence of voltages. 



magnetization using a self-consistent expression for the 
Keldysh Green's function. Parcollet and Hooley also ob- 
tain the finite field result M(B) = \h(B, V), which sug- 
gests that the voltage-dependent fluctuation-dissipation 
theorem is independent of field at weak coupling. 

Clearly, although this is beyond the scope of the cur- 
rent paper, the approach taken above can be extended 
to higher orders. An interesting question that this may 
help answer is whether the coherence of the Kondo effect 
is preserved at high voltage bias for the (physi cal) anti- 
ferromagnetically coupled quantum dot [HI fl?l ITsL ITflj . 

One of the enticing possibilities that this method offers 
is that of extension to more complex, multi-impurity or 
even lattice spin problems. The proportionality between 
spin and Majorana fermions can be extended to these 
cases, merely by introducing an independent Majorana 
fermion $j for each spin site, and writing ffj = 2<f>jSj. 
The generalization of (3) to a lattice is then 



identical result for J\\i = J was obtained in ref. [l4| . 
We can also read off the voltage-dependent fluctuation- 
dissipation relation (see Fig. 0), 



h(u) = - 



(1 + a)u 



[4>{oj + eV) + 4>{uj - eV)] + ot(j>{u) ' 



(36) 



where <f)(x) = x coth(x/2T), so that by (|17|l the static 
susceptibility is given by 



X(V,T) 



1 



(1 + a) 



AT [(eV/2T) coth(eV/2T) + a] 



(37) 



corresponding to a voltage-dependent Curie suscepti- 
bility. This result confirms earlier results of Parcollet 
and Hooley 0], obtained by a direct calculation of the 



i 



2Z, 



(38) 



where Zij — 2{&i$j). The quantity Z^ is a constant of 
motion that acts as a type of Z 2 gauge field. Closely 
related identities have recently been used to solve an 
anisotropic Heisenberg model on a honeycomb lattice 
[2(i|. The extension of these ideas to a Kondo lattice 
model, and its possible link to Z 2 gauge theories [2l| 
may be of particular interest to future research. 
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